N T I S ReP 
BY  PERM 
INFORMAll 


I 


^.-jcriON 

hI-»N  OF 

N CANADA 


Communications 
^ Research 
oo  Centre 


CO 

o 

< 

O 

<c 


CL. 

- -I 

I 

l 

O 

CD 

/ 

LU 
■ — " 1 

L 

c= 

2 

> uZ 

m 

■ 

AN  ITERATIVE  JECHNIQUE  FOR  L 

— —rssr 

DESIGNING  IINITEJIPULSE  RESPONSE 

JHEBYSHEV  IUIILTERS 

"WITH  NONLINEAR  DELAY  / 


by 


Ijq]  R.W./ Herring 


This  document  has  boon  a;-,  moved 
for  public  relcoc?  cr.d  scb;  it a 
distribution  is  unlimited. 


1 * 


% 


DEPARTMENT  OF  COMMUNICATIONS 
MINISTERE  DES  COMMUNICATIONS 


CRCLTECHNICAL  NOTE, (NO.  691 


l 


CAW 


>8  07  2$ 


MAR 


78 


.U 


Thli  work  wu  tpontored  by  the  Department  of  Netlonal  Defence, 
Reeeerch  end  Development  Brench  under  Project  No.  33C69. 


COMMUNICATIONS  RESEARCH  CENTRE 

DEPARTMENT  OF  COMMUNICATIONS 
CANADA 


AN  ITERATIVE  TECHNIQUE  FOR  DESIGNING  FINITE  IMPULSE  RESPONSE  CHEBYSHEV  MTI 
FILTERS  WITH  NONLINEAR  PHASE  DELAY 


CRC  TECHNICAL  NOTE  NO.  691  March  1978 

OTTAWA 


This  work  was  sponsored  by  the  Department  of  National  Defence,  Retearch  and  Development  Branch  under  Project  No.  33C69. 


TABLE  OF  CONTENTS 


ABSTRACT  1 

1.  INTRODUCTION  1 

2.  BENEFITS  OF  NONLINEAR  PHASE  2 

3.  OUTLINE  OF  DESIGN  PROCEDURE  • 9 

4.  DESIGNING  THE  PROTOTYPE  FILTER  10 

5.  GAIN  NORMALIZATION 13 

6.  COMPUTER  PROGRAMS  13 

7.  SUMMARY 14 

8.  ACKNOWLEDGEMENT 14 

9.  REFERENCES 14 

APPENDIX  A - PROGRAM  MPMTIPSF  16 

APPENDIX  B - PROGRAM  MPMTIPSFO  25 

APPENDIX  C - SUBROUTINES 34 


AN  ITERATIVE  TECHNIQUE  FOR  DESIGNING  FINITE  IMPULSE  RESPONSE  CHEBYSHEV  MTI 

FILTERS  WITH  NONLINEAR  PHASE  DELAY 


R.W.  Herring 


ABSTRACT 

-i 

A technique  for  designing  finite-impulse-response  (FIR)  equi ripple 
high-pass  digital  filters  with  nonlinear  phase,  suitable  for  use  in  radar 
moving-target-indication  (MTI)  systems,  is  described.  The 
optimization  of  such  filters  in  both  the  time  and  frequency  domains 
is  discussed,  and  comparisons  are  made  between  the  spectral 
per*  mance  of  nonlinear-phase  filters  and  that  of  I inear  phase  filters 
designed  to  meet  the  same  stopband  bandwidth  and  attenuation 
specifications.  In  92  out  of  the  94  cases  exam'ned,  the  transition 
bandwidth  between  the  stopband  is  narrower  fo.  me  nonlinear-phase 
filter.  In  the  MTI  application,  this  characteristic  makes  it  possible  to 
detect  targets  of  lower  radial  velocity.  Listings  of  Fortran  programs 
for  designing  these  nonlinear-phase  filters  are  included. 

F~  - ■ 

1.  INTRODUCTION 

The  object  of  this  note  is  to  describe  a technique  for  designing 
finite-impulse-response  (FIR)  digital  filters  with  nonlinear  phase-delays 
for  use  in  radar  moving-target-indication  (Mil)  systems.  The  optimization 
of  such  filters  in  both  the  time  and  frequency  domains  is  discussed,  and 
comparisons  are  made  between  the  spectral  performance  of  nonlinear-phase 
filters  and  that  of  linear-phase  filters  designed  to  meet  the  same  stopband 
bandwidth  and  attenuation  specifications.  In  92  out  of  94  cases  examined, 
the  transition  bandwidth  between  the  filter  stopband  anu  passband  is  narrower 
for  the  nonlinear-phase  filter.  Ir  the  MTI  application,  this  characteristic 
makes  it  possible  to  detect  targets  of  lower  radial  velocity. 

It  is  well  known  from  digital  signal  theory  (e.g.,  [1])  that  delaying 
a sampled  data  sequence  by  a fixed  amount  imposes  a shift  or  delay  in  the 
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relative  phases  of  its  spectral  components  which  is  linear  with  frequency. 

One  of  the  great  advantages  of  FIR  filters  is  that  they  can  be  designed  to 
have  such  a linear  phase-delay,  so  that  their  only  effects  are  to  modify  the 
magnitudes  of  the  spectral  components  and  to  delay  the  sequence  by  a fixed 
amount. 

In  contrast,  infinite-impulse-response  (HR)  filters  and  nonlinear-phase 
FIR  filters  have  nonlinear  phase  versus  frequency  characteristics.  The  use  of 
nonlinear-phase  filters  may  increase  the  complexity  of  any  coherent  post  MTI 
signal  processing,  but  since  the  phase  characteristics  of  these  filters  are 
deterministic,  any  undesirable  effects  due  to  the  phase  nonlinearities  can  be 
compensated . 

Houts  and  Burlage  [2]  have  described  the  advantages  to  be  derived  by 
using  Chebyshev  equiripple  FIR  filters  in  MTI  systems,  and  they  have  published 
computer  programs  [3]  for  the  design  and  evaluation  of  Chebyshev  filters 
having  linear  phase-delay.  Their  design  procedure  is  based  on  a computer 
program  for  designing  equiripple  FIR  linear  phase  digital  filters  [4]  using 
the  Remez  exchange  algorithm. 


2.  BENEFITS  OF  NONLINEAR  PHASE 

The  removal  of  the  linear-phase  constraint  can  result  in  improved  MTI 
performance  in  both  the  spectral  and  the  time  domains.  The  parameters  of 
interest  in  MTI  filter  design  and  defined  in  Table  1 and  depicted  in  Figure 
1.  The  standard  definitions  of  equiripple  FIR  filter  parameters  (e.g.,  [5]) 
are  defined  in  Table  2 and  depicted  in  Figure  2.  Note  that  the  standard 
definitions  refer  to  a filter  designed  to  have  unity  gain  in  the  passband, 
whereas  for  MTI  applications,  it  may  be  desirable  to  have  non-unity  passband 
gain  in  order  to  have  0 dB  white  noise  power  gain  and/or  0 d3  minimum  pass- 
band  gain. 

The  improved  spectral  performance  obtainable  from  nonlinear-phase 
equiripple  FIR  filters  relative  to  linear-phase  equiripple  FIR  filters  can 
be  realized  in  three  different  ways.  First,  smaller  ripples  in  either  the 
stopband  or  passband  ripples  (or  both)  can  be  achieved  for  given  values  of 
^STOP’  ^ASS  anc*  Smaller  passband  ripples  mean  increased  Agg  and  thus 
greater  clutter  suppression.  Decreased  Rpg  means  that  a higher  detection 
threshold  can  be  used  without  a loss  of  target  visibility  due  to  signal 
attenuation  in  the  filter  passband. 

Second,  the  transition  band,  or  the  band  of  frequencies  between  fgTOP 
and  fpASS>  can  narrower.  Such  a narrowing  of  the  transition  band  is 

useful  when  enhanced  low-velocity  target  visibility  is  desired  in  the  presence 
of  stationary  clutter  of  finite  bandwidth. 

Third,  the  performance  of  a given  linear-phase  filter  can  be  approxi- 
mated, except  for  phase,  by  a nonlinear-phase  filter  with  smaller  N.  Thus 
enhanced  incoherent  integration  gain  can  be  achieved  from  a given  fixed 
number  of  radar  pulses,  since  a greater  number  of  independent  filtered  output 
pulses  are  then  available  for  integration  [3] . 
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Definitions  of  Equiripple  MTI  Filter  Parameters 
N number  of  weights  in  filter  impulse  response 

h(n)  (0  < n < N-1)  filter  impulse  response 

f Ooppler  frequency  (Hz) 

|H(f)|  absolute  magnitude  of  the  MTI  filter  response 

H(f)  complex  magnitude  of  the  MTI  filter  response 

fgTQp  frequency  of  upper  edge  of  filter  stopband  (Hz) 

fpAss  frequency  of  lower  edge  of  filter  passband  (Hz) 

fpR  radar  pulse-repetition  frequency  (Hz) 

Rpg  peak-to-peak  gain  ripple  in  filter  passband  (dB) 

Agg  minimum  filter  attenuation  - measured  from  bottoms  of  passband  ripples  to  peaks  of  stopband 

ripples  (dB) 

GyVNP  white-noise-power  gain  (dB) 

GpBM  minimum  gain  in  filter  passband  (dB) 


|h (f )|  (dB) 


Figure  1.  Definition  of  Equiripple  MTI  Filter  Parameters 


TABLE  2 


amplitude  of  passband  ripple  (linear) 
C2  amplitude  of  stopband  ripple  (linear) 
(other  definitions  as  in  Table  1) 


|H(«| 


Figure  2.  Standard  Definition  of  Equiripple  High-Pass  Filter  Parameters 

In  Tables  3 and  4 there  can  be  found  examples  showing  each  of  the  three 
forms  of  benefit  described  above.  The  filters  summarized  in  Table  3 were 
designed  to  have  AgB  = 50  dB  and  values  of  fsxOP  ranging  from  50  Hz  to  400  Hz. 
Those  summarized  in  Table  4 were  designed  to  have  ASB  = 30  dB  and  values  of 
fSTOp  ranging  from  100  Hz  to  400  Hz.  In  both  tables  fpB  = 2500  Hz.  For 
some  sets  of  design  parameters  no  results  are  shown,  since  such  filters 
could  not  be  designed  subject  to  the  constraints  Gwnp  = 0 dB  and  GpBM  = 0 dB 
( vide  Tables  3a  and  3e) . 

The  relaxation  of  the  linear- phase  constraint  also  allows  some  latitude 
in  the  choice  of  the  filter  impulse-response  function  (IRF) , because  there 
are  several  possible  IRFs  corresponding  to  a particular  spectral  amplitude 
function.  Each  of  these  IRFs  corresponds  to  a different  phase-delay  charac- 
teristic, but  since  phase  delay  is  of  no  concern,  it  now  becomes  possible  to 
select  a particular  IRF  on  the  basis  of  its  time  domain  characteristics. 

Three  possible  criteria  for  this  selection  are:  (I)  a mini-max  criterion, 
which  attempts  to  equalize  the  magnitudes  of  the  IRF  weights  by  selecting 
that  IRF  having  the  smallest  value  of  the  ratio  of  the  largest  to  smallest 
weights;  (II)  a criterion  which  minimizes  susceptibility  to  numerical  over- 
flow by  selecting  that  IRF  having  the  minimum  absolute  sum  of  its  weights, 
or  (III)  another  criterion  which  attempts  to  equalize  the  magnitude  of  the 
IRF  weight  by  selecting  that  IRF  having  the  minimum  variance  of  the  magni- 
tudes of  its  weights.  Criteria  I and  III  should  therefore  select  IRFs  which 
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are  less  susceptible  to  disruption  by  impulsive  noise,  whereas  Criterion  II 
is  based  on  minimizing  susceptibility  to  numerical  overflow.  The  evaluation 
of  the  suitability  of  these  criteria  or  the  proposal  of  others  remains  as  a 
topic  requiring  further  investigation. 


TABLE  3a 

Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE 


NONLINEAR  PHASE 


mber  of 
lights  N 

Passband  Low-. 
Edge  (Hz)  fpASS 

Passband  Ripple 
<dB)  Rp0 

Passband  Lower 
Edge  (Hz)fpASS 

Passband  Ripple 
<dB>  RpB 

Reduction  of 
Passband  Lower 
Edge  (Hz) 

Reduction 
in  Passband 
Ripple  (dB) 

5 

477 

2 71 

_ 

_ 

„ 

__ 

6 

537 

3.30 

390 

2.12 

147 

1.18 

7 

392 

2.48 

376 

2.07 

16 

041 

8 

421 

2.52 

347 

1.93 

74 

0.59 

9 

386 

2.16 

319 

1.79 

67 

0.37 

10 

350 

2.05 

292 

1.64 

58 

0.41 

11 

352 

2.00 

270 

1.51 

82 

0.49 

12 

301 

1.73 

251 

1.40 

50 

0.33 

13 

317 

1.82 

234 

1.30 

83 

0.52 

14 

265 

1.51 

219 

1.22 

46 

0.29 

15 

286 

1.65 

206 

1.14 

80 

0.51 

16 

237 

1.34 

195 

1.08 

42 

0.26 

'stop  " 

25  Hz  fpR 

= 2500  Hz 

ASB  -50dB 

GWNP  = 

0dB  gpbm 

= 0dB 

TABLE  3b 

Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE 


NONLINEAR  PHASE 


Number  of 
Weights  N 

Passband  Lower 
Edge  (Hz)  fpASS 

Passband  Ripple 
(dB)  RpB 

Passband  Lower 
Edge  (Hz)  fpASS 

Passband  Ripple 
(dB)  Rpe 

Reduction  of 
Passband  Lower 
Edge  (Hz) 

Reduction 
in  Passband 
Ripple  (dB) 

5 

718 

4.30 

576 

3.48 

142 

0.82 

6 

539 

3.33 

482 

2.87 

57 

C.46 

7 

556 

3.44 

416 

2.45 

140 

0.99 

8 

423 

2.54 

367 

2.14 

56 

0.40 

9 

455 

2.80 

332 

1.90 

123 

0.90 

10 

352 

2.07 

300 

1.72 

52 

0.35 

11 

386 

2.35 

282 

1.60 

104 

0.75 

12 

303 

1.76 

280 

1.60 

23 

0.16 

13 

337 

2.03 

274 

1.57 

63 

0.46 

a 

14 

288 

1.66 

263 

1.52 

25 

0.14 

- 

15 

300 

1.79 

252 

1.46 

48 

0.33 

16 

284 

1.65 

240 

1.39 

44 

0.26 

‘stop 

50  Hz  <PR  = 

2500  Hz 

Agg  = 50  dB 

gwpn  ~ 0 dB 

gpbm 

— 0 dB 

<Mj* « , w*vvr~ < ^ ■ 
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TABLE  3c 

Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE  NONLINEAR  PHASE 


Reduction  of  Reduction 

Number  of  Passbend  Lower  Passband  Ripple  Passband  Lower  Passband  Ripple  Passband  Lower  in  Passband 


Weights  N 

Edge  (Hz)  fpASS 

<dB)  RpB 

Edge  (Hz)  fpASS 

<dB)  RpB 

Edge  (Hz) 

Ripple  (dB) 

5 

793 

5.38 

591 

3.67 

202 

1.71 

6 

542 

3.37 

532 

3.37 

10 

0.00 

7 

583 

3.78 

511 

3.13 

72 

0.65 

8 

537 

3.33 

462 

2.84 

75 

0.49 

9 

469 

2.96 

417 

2.56 

62 

0.40 

10 

477 

3.01 

380 

2.31 

97 

0.70 

11 

397 

2.47 

349 

2.11 

48 

0 36 

12 

419 

2.64 

332 

2.01 

87 

0.53 

13 

346 

2.13 

330 

2.00 

16 

0.13 

14 

372 

2.33 

321 

1.95 

51 

0.28 

15 

340 

2.09 

308 

1.88 

32 

0.21 

16 

335 

2.09 

293 

1.79 

42 

0.30 

’stop  = 

100  Hz  fpR 

= 2500  Hz 

ASB  = 50dB 

GWNP=0dB 

GPBM=0dB 

I 

!‘ 

I 

t; 

i; 


TABLE  3d 

Comparative  Examples  of  Equiripple  Finitc-Impulse-Response  MT!  Filter  Characteristics 


LINEAR  PHASE 


NONLINEAR  PHASE 


mber  of 

Passband  Lower 

Passband  Ripple 

Passband  Lower 

Passband  Ripple 

Reduction  of 

Passband  Lower 

Reduction 

in  Passband 

lights  N 

Edge  (Hz)  fpASS 

(dB)  RpB 

Edge  (Hz)  fpASS 

(dB)  Rp0 

Edge  (Hz) 

Ripple  (dB) 

5 

810 

5.65 

807 

5.61 

3 

0.04 

6 

818 

5.76 

648 

4.28 

170 

1.40 

7 

654 

3.55 

634 

4.16 

20 

-0.61 

8 

627 

4.26 

580 

3.81 

47 

0.45 

9 

617 

4.14 

522 

3.40 

95 

0.  »4 

10 

519 

3.45 

479 

3.30 

40 

0.15 

11 

536 

3.60 

479 

3.10 

57 

0.50 

12 

493 

3.25 

458 

2.97 

35 

0.28 

13 

471 

3.14 

431 

2.80 

40 

0.34 

14 

472 

3.14 

407 

2.64 

65 

0 50 

15 

423 

2.80 

405 

2.62 

18 

0.18 

16 

436 

2.90 

397 

2.58 

39 

0.32 

’stop 

200  Hz  fpR  : 

- 2500  Hz 

AS8  =50dB 

gwnp" 

0 dB  Gj 

PBM  ~ 0 dB 

j 
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MSif  3e 

Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE 

NONLINEAR  PHASE 

Number  of 

Pass  band  Lower 

Passband  Ripple 

Passband  Lower 

Passband  Ripple 

Reduction  of 
Passband  Lower 

Reduction 
in  Passband 

Weights  N 

Edge  (Hz)  *pASS 

(dB)  RpB 

Edge  (Hz)  fpA$s 

(dB)  RpB 

Edge  (Hz) 

Ripple  (dB) 

b 

- 

- 

- 

- 

- 

- 

6 

- 

- 

- 

- 

- 

1 

- 

- 

871 

6.66 

- 

- 

8 

873 

6.65 

789 

4.89 

84 

1.76 

y 

815 

5.58 

754 

5.50 

61 

0.08 

10 

716 

5.28 

684 

5.58 

32 

-0.30 

11 

723 

5.35 

686 

4.92 

37 

0.43 

12 

703 

5.16 

646 

4.61 

57 

0.55 

13 

636 

5.39 

636 

4.52 

0 

0.87 

14 

647 

4.74 

618 

4.39 

29 

0.35 

15 

640 

4.68 

596 

4.07 

44 

0.61 

16 

607 

3.97 

593 

4.20 

14 

-0.23 

'stop  = 

400  Hz  fpR 

= 2500  Hz 

ASB  = 50dB 

gwnp  = 

OdB  GpBM=0dB 

TABLE  4a 

Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE 


NONLINEAR  PHASE 


Number  of 
Weights  N 

Passband  Lower 
Edge  (Hzl  fpASS 

Passband  Ripple 
(dB)  RpB 

Passband  Lower 
Edge  (Hz)  fpASS 

Passband  Ripple 
(dB)  RpB 

Reduction  of 
Passband  Lower 
Edge  (Hz) 

Reduction 
in  Passband 
Ripple  (dB) 

5 

492 

2.65 

489 

2.78 

3 

-0.13 

6 

534 

3.27 

455 

2.61 

79 

0.66 

7 

485 

2.79 

407 

2.35 

78 

0.44 

8 

422 

2.53 

365 

2.12 

57 

041 

9 

428 

2.53 

330 

1.91 

98 

0.62 

10 

353 

2.09 

302 

1.75 

51 

0.34 

11 

374 

2.23 

279 

1.62 

95 

0.61 

12 

306 

1.80 

270 

1.56 

36 

0.24 

13 

332 

1.98 

268 

1.55 

64 

0.43 

14 

278 

1.49 

262 

1.53 

16 

-0.04 

15 

298 

1.78 

253 

1.48 

45 

0.30 

16 

275 

1.61 

243 

1.43 

32 

0.18 

'stop  ~ 

100  Hz  fpR 

= 2500  Hz 

Asb  = 30  dB 

°WNP  ” 

0<1B  gpbm 

= 0dB 
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TABLE  4b 


Comparative  Examples  of  Equiripple  Finite-Impulsc-Response  MTI  Filter  Characteristics 

LINEAR  PHASE  NONLINEAR  PHASE 


Number  of 


1 1 

Weights  N 

Edge  (Hz)  fpAS$ 

(dB)  Rps 

1 

5 

766 

4.97 

II  , 

6 

549 

3.46 

f|  i 

7 

479 

3.73 

fL 

8 

521 

3.28 

I 1 

9 

474 

3.03 

W 

10 

479 

3.05 

1' 

11 

409 

2.61 

| 

12 

428 

2.74 

j 

13 

395 

2.52 

14 

386 

2.48 

p/ 
pf  . 

15 

384 

2.46 

16 

353 

2.28 

Passband  Lower 

Passband  Ripple 

Reduction  of 
Passband  Lower 

Reduction 
in  Passband 

Edge(Hz)fpASS 

fdSIRpB 

Edge  (Hz) 

Ripple  (dB) 

594 

3.71 

172 

1.26 

503 

3.64 

46 

-0.18 

505 

3.12 

74 

0.61 

469 

2.91 

53 

0.37 

428 

2.67 

46 

0.36 

394 

2.47 

85 

0.58 

387 

2.42 

22 

0.19 

380 

2.39 

48 

0.35 

366 

2.31 

29 

0.21 

349 

2.21 

37 

0.27 

334 

2.12 

50 

0.34 

331 

2.11 

22 

0.17 

fSTOP  ~ 200  Hz 


fpn  = 2500  Hz 


Agg ’ 30  dB 


GWNP=0dB 


gpbm  0 dB 


Number  of 
Weights  N 


TABLE  4c 


Comparative  Examples  of  Equiripple  Finite-Impulse-Response  MTI  Filter  Characteristics 


LINEAR  PHASE 


Pass  band  Lower 
EdgelH*)fpASS 


Passbend  Ripple 
(dB)  R_„ 


NON  LINEAR  PHASE 
nd  Lower  Passbaud  Ripple 


Passband  Lower 
Edge(Hz)fpASS 


Reduction  of 
Passband  Lower 


Reduction 


(dB>  Rpb 

Edge  (Hz) 

Ripple  1 

5.84 

0 

0.01 

4.66 

70 

1.29 

4.85 

76 

0.21 

4.97 

-30 

-0.31 

4.35 

39 

0.40 

4.10 

52 

0.45 

4.00 

9 

0.17 

3.92 

30 

0.32 

3.77 

41 

0.38 

3.75 

6 

0.16 

7.68 

25 

0.28 

3.55 

30 

0.35 

fST0P ~ 400  H? 


(pR  = 2500  Hz 


Agg  ^ 30  dB 


^WN ? ^ 


GPBM  = 0d8 


m *v 
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3.  OUTLINE  OF  THE  DESIGN  PROCEDURE 

The  technique  for  designing  nonlinear-phase  FIR  filters  is  based  on  a 
procedure  first  suggested  in  [6].  This  procedure  involves  the  design  of  a 
prototype  linear-phase  FIR  filter  having  an  IRF  hp(n)  of  length  (2N-I)  and  an 
amplitude  response  in  the  frequency  domain  Hp(f)  equal  to  the  square  of  the 
magnitude  of  the  desired  frequency  response.  It  is  shown  below  that  this 
prototype  filter  cannot  be  designed  directly  from  the  desired  parameters  of 
the  nonlinear-phase  filter,  but  that  an  iterative  technique  must  be  used. 

A property  of  linear-phase  FIR  filters  with  real-valued  IRFs  is  that  if 
z0  is  a zero  of  the  IRF,  then  so  are  z0-f,  z*  and  (z0--*-)*  where  * denotes 
complex  conjugate  ([1];  page  159).  Note  that  if  zQ  is  real,  then  zQ  = z*  and 
also,  if  jz0j  = 1,  then  zQ~^  = z*,  so  that  roots  can  be  real  and  single,  or 
occur  in  conjugate  pairs  (if  |zQ|  = 1)  or  in  reciprocal  pairs  (if  zQ  is  real), 
or  in  conjugate-reciprocal  quartets  (if  |z0|  ^ 1 and  zQ  complex).  In  general, 
of  the  2(N-1)  zeroes  of  the  (2N-l)-length  prototype  filter,  there  are  k pairs 
of  reciprocal  real  zeroes,  £ quartets  of  conjugate  reciprocal  zeroes  and 
m = (N-k-2£-l)  pairs  of  double  zeroes  on  the  unit  circle.  To  extract  an  IRF 
h(n)  of  length  N,  it  is  necessary  to  discard  one  zero  of  each  of  the  m pairs 
of  double  zeroes,  one  zero  of  each  of  the  k pairs  of  reciprocal  zeroes,  and 
one  pair  of  conjugate  zeroes  from  each  of  the  £ quartets  of  conjugate- 
reciprocal  zeroes.  This  leaves  a set  of  (N-l)  zeroes  which  is  expanded  to 
produce  a real-valued  IRF  of  length  N. 

A set  of  M IRFs  of  length  N can  thus  be  derived  from  the  (2N-l)-length 
prototype  filter,  where 

M = 2^C+*'-^ 

Note  that  one  zero  or  pair  of  conjugate  zeroes  can  be  arbitrarily  chosen  to 
lie  inside  or  outside  the  unit  circle,  since  the  only  effect  of  this  choice 
is  to  reverse  the  IRFs  in  the  time  domain;  i.e.,  h(n)  is  replaced  by  h(N-l-n) 
for  0 <_  n _<  N-l . 

An  upper  bound  on  M as  a function  of  N is  given  by 

M _ ,([(N-2)/2]-l) 

where  [x]  denotes  the  largest  integer  less  than  x,  so  it  can  be  seen  that  the 
set  of  IRFs  to  be  examined  can  contain  of  order  210  for  N=24.  Hence  the 
optimization  techniques  described  above  can  become  quite  time-consuming  for 
filters  of  such  length.  This  limitation,  however,  should  not  preclude  the 
use  of  these  procedures  for  designing  filters  of  the  lengths  usually  consid- 
ered for  MTI  applications  (e.g.,  N <_  20;  so  that  M <_  256). 


4.  DESIGNING  THE  PROTOTYPE  FILTER 

In  order  to  make  use  of  the  standard  filter  design  algorithms,  it  is 
first  necessary  to  define  the  filter  ripple  parameters  (Rpg,Agg)  of  interest 
to  radar  MTI  designers  in  terms  of  the  standard  linear  ripple  parameters 


V A -A-  .'~i*  . 


C-yVV*-- «•*> 


uKLSsmiLM 


(e^.e^).  It  can  easily  be  shown  from  Figures  1 and  2 that 


lo(V20)  1 


cl  = (R  /20)  ~ 
10  PB  +1 


1 - c. 


2 


In  the  standard  design  procedure  for  linear-phase  FIR  filters  [3],  if 


N,  fST0P’  fPASS  and  fPR  are  sPecified»  then  only  the  ratio 


W = C2./£2 


can  be  specified  as  a free  parameter.  This  restriction  can  be  circumvented 
by  allowing  fp^ss  t0  varied  in  an  iterative  manner  until  that  value  for 
fpAss  Is  found  which  gives  the  desired  values  for  S£  and  hence  e-^  [3], 


A similar  iterative  procedure  is  necessary  in  the  design  of  nonlinear- 
phase  FIR  filters.  For  the  design  of  the  linear-phase  prototype  filter  H0(f) 


(see  Figures  3(a)-3(c)),  it  is  necessary  to  specify  2N-1,  ^STOP’  fpR>  ^1  and 


6 2,  where  6i  and  $2  have  to  be  specified  in  terms  of  and  C2-  The  outline 
of  the  scheme  for  relating  the  6’s  and  e's  is  pictured  in  Figures  3(a)-3(c) 
and  is  similar  to  that  of  [6]  except  for  one  detail  pointed  out  below.  For 
clarity  the  scheme  is  described  progressing  from  the  original  prototype  filter 
H0(f)  (Figure  3(a))  to  the  intermediate  filter  Hr(f)  (Figure  3(b))  to  the 
final  prototype  filter  Hp(f)  = | H ( f ) | (Figure  3(c)),  In  practice,  the  actual 
progression  is  from  specifying  Hp(f)  to  specifying  H0(f)  in  terms  of  Hp(f), 
since  H0(f)  is  the  filter  which  is  actually  designed  using  the  algorithm 
of  [4] . 


H (f)  Is  related  to  H (f)  by  the  transformation 

A C 


H (f ) = H (f)  + 6„. 
I o i 


In  the  time  domain  this  is  equivalent  to 


hx(n)  = 


h (n) , 


1 < n < N-l 


h (n)  + 5 , n = 0. 
o z 


H (f)  is  in  turn  related  to  H_(f)  by  the  transformation 
P I 


Hp(f)  = K Hx(f) 


which  becomes  in  the  time  domain 


(7) 


•1 


n)  = K hjCn), 


-N+l  < n < N-l 


or.  It  is  in  the  determination  of  the  factor  K that 
from  that  of  [6],  for  it  can  be  noted  that  a scaling 
es  not  produce  a filter  | H (f ) | with  equal  (linear) 
ns  above  and  below  unity  gain. 


ameters  Rpg  and  Agg  and  thus  ep  and  c£  which  are 
esignar,  it  is  necessary  to  invert  the  above 
o derive  expressions  for  6^  and  <$2  in  terms  of 
it  is  a straightforward  matter  to  use  the  linear- 


It  is  now  a simple  matter  to  derive  hp(n)  from  hQ(n)  in  terms  of  these  factors 
as  outlined  above,  and  then  to  derive  a nonlinear-phase  N-length  IRF  h(n). 


5.  GAIN  NORMALIZATION  .j 

It  is  often  useful  to  design  MTI  filters  which  have  0 dB  white-noise  | 

gain.  This  means  that  the  white-noise  power  in  the  filtered  signal  is  the  J 

same  as  that  in  the  unfiltered  input  signal,  so  that  in  the  absence  of  clutter,  j 

the  false-alarm  probability  Pp^  is  not  altered.  Such  normalization  is  easily  i 

accomplished  [3]  by  invoking  Parseval’s  theorem  -j 

n-i  r0.5  fp_  | 

g2  = 1 |h(n) | 2 = 2/  PR  | H(f ) | 2 df  (14) 

n~o  JQ  i 

to  compute  the  white  noise  power  gain  g2  by  means  of  a simple  summation.  . 

Scaling  h(n)  by  g--*-  produces  a filter  with  unity  white  noise  power  gain 
(^WNP  = 0 dB)  and  leaves  Acr  and  RpR  unaltered.  Note  that  this  procedure  is 
valid  for  any  FIR  filter.  ]j 

It  can  also  be  useful  to  design  a filter  which  has  both  G^p  = 0 dB  and  ii 

a minimum  gain  of  unity  in  its  passband  (GpgM  = 0 dB) . This  means  that  the  v 

probability  of  detection  Pp  in  the  passband  will  not  be  reduced  by  filter 
attenuation  in  these  passband  ripples.  This  goal  can  be  achieved  by  an  >j 

iterative  procedure.  In  this  case  it  is  the  passband  ripple  that  is  altered  j 

until  the  desired  characteristics  are  obtained.  The  procedure  is  to  specify  ’ 

N,  Agg  and  fgg,  select  an  arbitrary  reasonable  value  for  Rpg  and  design  a l| 

filter  using  the  iterative  procedure  described  above.  Based  on  whether  this 
filter  has  GpBM  greater  or  less  than  unity  gain,  the  value  for  Rpg  is  ? 

increased  or  decreased  respectively  and  another  iteration  is  carried  out. 

'This  procedure  is  repeated  until  convergence  to  Gpg^  = 0 dB  is  achieved  to  | 

within  an  acceptable  tolerance.  >: 

f 

3 

l\ 

6.  COMPUTER  PROGRAMS  I 


The  program  Ml’MTIPSF  generates  nonlinear-phase  FIR  MTI  filters  for 
||.  specified  values  of  N,  fgTOp>  fpR  Agg  and  Rp^.  Parameters  to  specify  the 


*?*«!'*'*?  ■ *-£i  ^<7^hV^.n  »'■'<*■'£&  rtv»  ■ ,-<"<*  V *•  cv^4*^/j  <^V* 
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grid  density  for  the  design  of  the  prototype  filter  [4]  and  to  select  the 
impulse-response  design  criterion  must  also  be  provided.  This  program  is 
listed  in  Appendix  A,  and  the  detailed  operating  instructions  are  given  there. 
Typical  times  required  to  design  a filter  using  a Xerox  Sigma  9 computer  range 
from  2.2  sec  for  N=5  to  28.4  sec  for  N-16. 

The  program  MPMTIPSFO  generates  nonlinear-phase  FIR  MTI  filters  with 
GpgM  = 0 dB.  The  same  parameters  as  for  the  program  MPfT'IPSF  must  be  specified 
by  the  user,  but  here  the  parameter  Rp^  is  used  only  as  an  initial  value  which 
is  then  modified  by  the  program  to  converge  to  Gpg^  = 0 dB.  Hence  the  choice 
of  a starting  value  for  Rpg  near  its  final  value  can  significantly  reduce  the 
time  required  to  design  a filter.  This  program  is  listed  in  Appendix  B. 

Typical  running  times  on  the  Sigma  9 can  be  1 to  10  or  more  times  as  long  as 
for  the  program  MPMTIPSF,  depending  on  how  good  an  estimate  for  the  initial 
value  of  Rv^  is  used. 

Appendix  C contains  the  listings  of  the  subroutines  required  by  the 
programs  MPMTIPSF  and  MPMTIPSFO.  The  subroutine  REMEZ  and  its  ancillary 
subroutines  have  not  been  included,  since  they  are  readily  accessible  else- 
wh  re  [4],  but  note  the  modifications  required  in  the  dimensions  of  the 
COMMON  block  variables. 


7.  SUMMARY 


It  has  been  shown  in  Tables  3 and  4 that,  for  the  typically  narrow  stop- 
bands  used  in  radar  MTI  filters,  nonlinear-phase  filters  can  be  designed  which 
offer  superior  visibility  for  low-velocity  targets,  relative  to  that  offered 
by  linear-phase  filters  designed  to  the  same  stopband  specifications.  It  has 
arso  been  pointed  out  that  the  time-domain  response  of  such  nonlinear-phase 
MTI  filters  can  be  optimized  without  altering  the  filter  power  response  in 
the  frequency  domain.  The  details  of  the  design  procedure  have  been  described, 
and  listings  of  Fortran  programs  for  implementing  the  procedure  have  been 
provided. 
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